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The hydrodynamic equations of superfluids for a weakly in- 
teracting Bose gas are generalized to include the effects of pe- 
riodic optical potentials produced by stationary laser beams. 
The new equations are characterized by a renormalized inter- 
action coupling constant and by an effective mass accounting 
for the inertia of the system along the laser direction. For 
large laser intensities the effective mass is directly related to 
the tunneling rate between two consecutive wells. The pre- 
dictions for the frequencies of the collective modes of a con- 
densate confined by a magnetic harmonic trap are discussed 
for both ID and 2D optical lattices and compared with recent 
experimental data. 

The experimental realization of optical lattices [Q-^ is 
stimulating new perpectives in the study of coherence 
phenomena in trapped Bose-Einstein condensates. A 
first direct measurement of the critical Josephson current 
has been recently obtained in ^ by studying the center 
of mass motion of a magnetically trapped gas in the pres- 
ence of a ID periodic optical potential. Under these con- 
ditions the propagation of collective modes is a genuine 
quantum effect produced by the tunneling through the 
barriers and by the superfluid behaviour associated with 
the coherence of the order parameter between different 
wells. The effect of the optical potential is to increase the 
inertia of the gas along the direction of the laser giving 
rise to a reduction of the frequency of the oscillation. 

The purpose of the present work is to investigate the 
collective oscillations of a magnetically trapped gas in the 
presence of ID and 2D optical lattices taking into account 
the effect of tunneling, the role of the mean field inter- 
action and the 3D nature of the sample. Under suitable 
conditions these effects can be described by properly gen- 
eralizing the hydrodynamic equations of superfiuids |^] . 

Let us assume that the gas, at T = 0, be trapped by an 
external potential given by the sum of a harmonic trap of 
magnetic origin Vho and of a stationary optical potential 
V^pt modulated along the z-axis. The resulting potential 
is given by 



where ujr 
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, ujz are the frequencies of the harmonic trap, 
q = 27r/A is fixed by the wavelength of the laser light cre- 
ating the stationary ID lattice wave, Er = fi^q^ /2m is 



the so called recoil energy and s is a dimensionless pa- 
rameter providing the intensity of the laser beam. The 
optical potential has periodicity d = 7r/g = A/2 along the 
z-axis. The case of a 2D lattice will be discussed later. In 
the following we will assume that the laser intensity be 
large enough to create many separated wells giving rise 
to an array of several condensates. Still, due to quan- 
tum tunneling, the overlap between the wave functions 
of two consecutive wells can be sufficient to ensure full 
coherence. In this case one is allowed to use the Gross- 
Pitaevskii (GP) theory for the order parameter to study 
both the equilibrium and the dynamic behaviour of the 
system at zero temperature Eventually, if the tun- 
nelling becomes too small, the fiuctuations of the relative 
phase between the condensates will destroy the coherence 
of the sample giving rise to new quantum configurations 
associated with the transition to a Mott insulator phase 
i|. 

In the presence of coherence it is natural to make the 
ansatz 



fc 



(2) 



for the order parameter in terms of a sum of many 
condensate wave-functions relative to each well. Here 
Sk{x,y) is the phase of the /c-component of the order 
parameter, while ^'fc and fk are real functions. We 
will make the further periodicity assumption fk{z) — 
fo{z — kd) where /o is localized at the origin. The above 
assumptions for ^ and fk are justified for relatively large 
values of s where the interwell barriers are significantly 
higher than the chemical potential. In this case the con- 
densate wave functions of different sites are well sepa- 
rated (tight binding approximation). 

Using the ansatz (|^) for the order parameter one finds 
the following result for the mean field expectation value of 
the effective Hamiltonian H = (p|/2to -t- Vexti^j)) + 
.9Ej<fc'5(i-j -rfc): 
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where in the two-body and in the magnetic interaction 
terms as well as in the radial kinetic energy we have 
ignored the overlap contributions arising from different 
wells. In the evaluation of the axial kinetic energy and of 
the optical potential term we have instead kept also the 
overlap terms originating from consecutive wells. These 
are proportional to the quantity 



5 = -2 dz—dJoiz)d^.foiz~d) + foiz)foiz-d)Vopt , 
J 2m 
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related to the tunneling rate and responsible for the oc- 
currence of Josephson effects. 

By setting Sk = (groundstate configuration), the 
variation of E with respect to /o yields the differential 
equation 



2m dz^ 



sEji sin qz 
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where Sq is introduced to ensure the normalization con- 
dition /^^^2 ^•^/o ~ 1 which implies that the functions 
^'fc are normalized to the number of atoms Nk occu- 
pying each site: J "^^dxdy = Nk- In eq. (|^) we have 
ignored the contribution arising from the two-body in- 
teraction. Estimates of show that this is a good ap- 
proximation already at moderately large s. We have also 
neglected the external magnetic potential which is justi- 
fied if ^/sEji » hux, hujy. Since in the following we are 
interested in the low energy excitations of the system we 
will always keep the function /o equal to the groundstate 
solution of (H). 

In order to discuss the macroscopic properties of the 
system, including its low energy dynamics, it is conve- 
nient to transform the discretized formalism described 
above into the one of continuum variables. This is ob- 
tained through the replacement — *■ (l/d) Jdz in the 
various terms of the energy. Through such a procedure 
one naturally introduces a smoothed or "macroscopic" 
density defined by 



nM{x,y,z) = 



(6) 



with z ~dk, and a smoothed phase S (Sk S{x,y,z)). 

By applying the smoothing procedure to eq. we ob- 
tain the following macroscopic expression for the energy 
functional 
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where we have introduced the renormalized coupling con- 
stant g — gd J fgdz , we have neglected quantum pressure 
terms originating from the radial term in the kinetic en- 
ergy and we have set '^k'i'k+i ^ = duM- We have 
also omitted some constant terms (first two terms in eq. 
(^)) which do not depend on ni\j or on S. 

With respect to the functional characterizing a trapped 
Bose gas in the absence of optical confinement, one no- 
tices two important differences: first the interaction cou- 
pling constant is renormalized due to the presence of the 
optical lattice. This is the result of the local compres- 
sion of the gas produced by the tight optical confinement 
which increases the repulsive effect of the interactions. 
Second the kinetic energy term along the z-direction has 
no longer the classical quadratic form as in the radial 
direction, but exhibits a periodic dependence on the gra- 
dient of the phase. By expanding this term for small 
gradients, which is the case in the study of small am- 
plitude oscillations, one derives a quadratic term of the 
form (ti^ /2m*) J dVuMidzS)"^ characterized by the effec- 
tive mass 

m mSd^ (5 TT^ 

(8) 



Notice that within the 



where S is defined by eq. ([_ 
employed approximation the value of S, and hence of m*, 
does not depend on the number of atoms, nor on the 
mean field interaction. 

The equilibrium density profile, obtained by minimiz- 
ing eq. with 5 = has the typical form of an inverted 
parabola 
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which conserves the aspect ratio of the original mag- 
netic trapping. The size of the condensate has instead 
increased since g > g. For large s the increase of the cou- 
pling constant can be large {g ^ s^/^ |^). However, since 
the radius of the sample scales like the 1/5-th power of 
g the resulting increase in the size of the system is not 
very spectacular (for s = 15 we find an increase of the 
size by ^ 20% for the experimental setting of [||). 

The functional (0) can be used to carry out dy- 
namic calculations. In this case one needs the action 
A = J dt [{H) — , with the second term given by 

ih{{d/dt)) — — J dVTinMS. The resulting equations of 
motion are obtained by imposing the stationarity condi- 
tion on the action with respect to arbitrary variations of 
the density n^/ and of the phase S. The equations take 
the form 

Ti Sd 
riM + —dr^ {nndr^S) -t- —d^ (riMsin [dd^S]) = , (10) 
m h 

nS + guM + V^ag + ^ {dr^Sf ~ Scos [dd,S] = . (11) 
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In particular, at equilibrium these equations reproduce 
result (|^) for the equilibrium density. Furthermore, 
Josephson-type oscillations are among those captured by 
eqs. (|o|) and To see this consider the case of a 

uniform gradient of the phase along z, dzS = Pz{t)/h, 
where Pz is a time-dependent parameter. From eqs. ( p^ ) 
and ( pT| ) one can then derive equations of motion for the 
center of mass Z{t) = J dVznm{t)/N and for the conju- 
gate momentum variable Pz P,pT[| 



fiZ — 5dsm 



= 0, 



Pz + mujiZ = 0, 



(12) 
(13) 



which have the typical Josephson form. 

In the limit of small oscillations the solutions of eqs. 
( p^ ) and ( pT| ) have the form n = ri°j -I- Sn{r)e'"^* with dn 
obeying the hydrodynamic equations: 



dr , Sn 



Pl. - Vho 



dzSn 



, (14) 



where /i = 3"m(0) is the chemical potential of the sample 
and n'^,j{0) is the equilibrium density (^ evaluated at 
the center. The solutions of (^4|) provide the low energy 
excitations of the system. In the absence of magnetic 
trapping one finds phonons propagating at the velocity 
c = \/.9^A/f/"^*i agreement with the result obtained 
in |13] for a ID array of Josephson junctions. In the 
presence of harmonic trapping the discretized frequencies 
of the time-dependent solutions of ( |l4|) do not depend 
on the value of the coupling constant. By applying the 
transformation z \/ m* /mz, one actually finds that 
the new frequencies are simply obtained from the results 
of 10] by replacing 



(15) 



For an elongated trap (lo^ — ojy — lo±_ ^ lJz) the 
lowest solutions are given by the center-of-mass mo- 
tion lub = -y/ m/m*LOz and by the quadrupole mode 
loq = y/5/2y/m/m*LUz. The center-of-mass frequency 
coincides with the value obtained from eqs. (|lj) and 
( |l3| ) in the limit of small oscillations. Concerning the 
quadrupole frequency we note that the occurrence of the 
factor ■y/5/2 is a non-trivial consequence of the mean 
field interaction predicted by the hydrodynamic theory 
of superfluids in the presence of harmonic trapping [Q. 
In addition to the low-lying axial motion the system ex- 
hibits radial oscillations at high frequency, of the order of 
uj± . The most important ones are the transverse breath- 
ing and quadrupole oscillations occuring at uj — 2lu± and 
UJ = ^/2uj± respectively. For elongated traps the frequen- 
cies of these modes should not be affected by the presence 
of the optical potential. Different scenarios are obtained 
for disc-shaped traps {wz » ^±)- The above results ap- 
ply to the linear regime of small oscillations. Eqs. (|l3) 



and (13) show that in the case of center-of-mass oscilla- 
tions, the linearity condition is achieved for initial dis- 
placements Aa; of the trap satisfying Ax « ^2S/mu!'^, 
a condition that becomes more and more severe as the 
laser intensity increases. For larger initial displacements 
the oscillation is described by the pendulum equations. 
For very large amplitudes the motion is however dynam- 
ically unstable |ll|,0. 

From the previous discussion it emerges that the ef- 
fective mass is the crucial parameter needed to pre- 
dict the value of the small amplitude collective frequen- 
cies. An estimate of m/m* can be made by neglecting 
the magnetic trapping as well as the role of the mean 
field interaction. Within this approximation the effec- 
tive mass is easily obtained from the excitation spec- 
trum of the Schrodinger equation for the ID Hamilto- 
man H = -{h'^/2m)d^/dz^ + sEr sir? qz, avoiding the 
explicit determination of the tunneling parameter (jj). 
One looks for solutions of the form e^P^^^ fp{z) where p is 
the quasi- momentum of the atom and fp{z) is a periodic 
function of period d. The resulting dispersion law e(p) 
provides, for small p, the effective mass according to the 
identification e{p) ~ Eq +p'^/2m*. The value of m/m*, 
which turns out to be a universal function of the inten- 
sity parameter s, has been evaluated for a wide range of 
values of s (see fig.||). These results for m* can be used 
to estimate the actual value of the collective frequencies. 
The method described here to calculate m* is expected 
to be reliable not only for very large laser intensities s 
when the tight binding approximation applies and the 
effective mass can be expressed in terms of the tunneling 
rate (see eqs. (|[),(^)), but also for smaller values of s. Of 
course for very small laser intensities, as in the experi- 
ment p^ ] , the determination of m* requires the inclusion 
of the mean field interaction and of the magnetic trap- 
ping through the explicit solution of the GP-equation. 

In fig. 1^ we compare our predictions for the frequencies 
of the center-of-mass motion with the recent experimen- 
tal data obtained in The comparison reveals good 
agreement with the experiments. Our results also agree 
well with those obtained from the numerical solution of 
the time-dependent GP-equation ||Tl| , [l^ . 

The above formalism is naturally generalized to in- 
clude a 2D optical lattice where the optical potential is 
Vopt = sEjisin^qx + sEjisin^qy. The actual potential 
now generates an array of ID condensates which has al- 
ready been the object of experimental studies [Q. For a 
2D-lattice the ansatz for the order parameter is [Ea 



$(r) 



*;c.,;c„(^)/fe.,fe„(x,2/)e'^'===-'=«(^) 



(16) 



In the TF-limit the groundstate smoothed density um = 
^fe fe /"^^ ^lil^ ^^'^ familiar form n\,j = (/i — Vho) /g 
with the redefined coupling constant g = g {dj dx f^)^ , 
where /o is still given by the solution of eq. (|[) and we 
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have used the same approximations as in the ID case. 

Also with regard to dynamics, one can proceed as for 
the ID lattice. One finds that the equations of motions, 
after linearization, take the form 



5n — dz 



Vy 



ho 



dzSn 



dr 



ho 



dr , Sn 



(17) 



The frequencies of the low energy collective modes are 
then obtained from those in the absence of the lattice 
[0 by simply replacing ujx y^mJrnFuJx and tOy — > 
■y/ m/m*ujy. For large laser intensities the value of m* 
coincides with the one calculated for the ID array. If 
ujz » ym/m* , ujy y^m/m*, the lowest energy solu- 
tions involve the motion in the x ^ y plane. The oscil- 
lations in the z-direction are instead fixed by the value 
oiujz- These include the center-of-mass motion {lu — luz) 
and the lowest compression mode {lu = ^/3w^) Qg. The 
frequency uj — \/3uJz coincides with the value obtained by 
directly applying the hydrodynamic theory to ID systems 
p6| , p7[ and reveals the ID nature of the tubes generated 
by the 2D lattice. If the radial trapping generated by the 
lattice becomes too strong the motion along the tubes 
can no longer be described by the mean field equations 
and one jumps into more correlated ID regimes [p^ . 
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FIG. 1. Effective mass as a function of the laser intensity 
s (see eq.(|l|)) calculated neglecting the effects of interaction 
and harmonic trapping. 
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FIG. 2. Frequency of the center-of-mass motion for a con- 
densate trapped by the combined magnetic and optical po- 
tential (^ as a function of the laser intensity. The circles and 
triangles are, respectively, the experimental and theoretical 
data of [^. The triangles have been obtained by evaluating 
the tunneling rate within a Gaussian approximation for the 
order parameter in each well The solid line refers to our 
theoretical prediction. 
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